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Abstract 

Three-dimensional viscous flow computations are 
presented for the F/A-18 forebody-LEX geometry. 
Solutions are obtained from an algorithm for the 
compressible Navier-Stokes equations which incorpo- 
rates an upwind-biased, flux-difference-splitting ap- 
proach along with longitudinally -patched grids. Re- 
sults are presented for both laminar and fully turbulent 
flow assumptions and include correlations with wind 
tunnel as well as flight-test results. A good quantita- 
tive agreement for the forebody surface pressure distri- 
bution is achieved between the turbulent computations 
and wind tunnel measurements at = 0,6. The 
computed turbulent surface flow patterns on the fore- 
body qualitatively agree well with in-flight surface flow 
patterns obtained on an F/A-18 aircraft at M 0 0 — 0,34* 
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Nomenclature 

drag coefficient, Drag/ qooSref 

lift coefficient, Lift/ q^Sref 

pressure coefficient, (p — Poo)/<Zod 

specific heat at constant pressure 

specific heat at constant volume 

mean aerodynamic chord of reference wing 

planform, 8.294 in, 

total energy per unit volume 

flux vectors 

total enthalpy per unit volume 

Jacobian of the coordinate transformation 

conductivity 

Mach number 

Prandtl number 

state vector, J -1 [p, pu, pv, pw> e\ T 
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q total velocity 

qco freestream dynamic pressure 

Rc Reynolds number based on ? 

S re j area of reference wing planform, 

extended to configuration centerline 
u, v, w body-axis Cartesian vel ocity components 
v* wall-friction velocity, \/r w /p 

yfs fraction of LEX exposed semispan 

inner-law variable, yv* jv 
a angle of attack, degrees 

7 ratio of specific heats, c p /c v 

fi viscosity 

v kinematic viscosity, pt.jp 

0 azimuthal angle 

17, f body-fitted coordinates 
p density 

r w wall shear stress 

Subscripts 

1 laminar 

t turbulent 

v viscous 

oo freestream reference conditions 


Introduction 

Advances in numerical solution methodology along 
with increased computer speed and capacity have made 
it feasible to seek numerical solutions to the three- 
dimensional Navier-Stokes equations for relatively com- 
plicated geometries. Prior computations for isolated 
wing 1,2 * 3 or body 4,5,6 geometries have demonstrated 
that accurate Navier-Stokes solutions are not limited 
to benign flow conditions; converged results can be 
obtained which include complex flow structures (e.g,, 
three-dimensional separation, shocks, vortices, etc.) as 
well as interactions among these structures. 

Applications to aircraft configurations have been 
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less forthcoming primarily due to complexities associ- 
ated with grid generation and surface definition along 
with the relatively large computational resources re- 
quired for these cases, None the less, recent progress 
has been shown for several cases including supersonic 
inviscid flow about the F/A-18 forebody-LBX and the 
SR-71 7 as well as transonic viscous flow about the F- 
16A S , Most recently, viscous flow computations about 
the ascent configuration of the space shuttle 9 have 
demonstrated good correlation with flight data at sub- 
sonic, transonic, and supersonic speeds. 

The present investigation is directed toward ap- 
plying an extended version of an implicit Navier-Stokes 
algorithm 1 ’ 3,6 to the F/A-18 forebody-LEX geometry 
at conditions suitable to the formation of separated 
subsonic and, in some cases, transonic flows. The ex- 
tension to the algorithm permits longitudinally-blocked 
grids which are necessary for accurately modeling the 
subject configuration. The F/A-18 has been selected 
for this study primarily due to the availability of current 
wind-tunnel 10 a3 well as flight research 11 data related 
to an ongoing high angle-of-attack research program 
being conducted by NASA. 

The computations are focused on two specific flow 
conditions in accordance with the recent wind-tunnel 
and flight-test research. The wind tunnel conditions 
are ~ 0.6, = 0.8 x 10 6 , and a — 20° which cor- 

respond to recent tests conducted at the David Taylor 
Research Center (DTRC) 7x10 foot tunnel with a 6% 
F/A-18 model. At these conditions the forebody flow 
could be transitional, hence both laminar as well as tur- 
bulent solutions are obtained. The flight conditions are 

~ 0.34, f? 5 = 13.5 X 10 6 , and a = 19° which corre- 
spond to recent flight tests of the NASA F/A-18 High 
Alpha Research Vehicle (HARV) conducted at the Dry- 
den Flight Research Facility. At these conditions only 
turbulent solutions are obtained, but grid sensitivity is 
addressed. 

Governing Equations 

The governing equations as well as computational 
method for the present investigation have been pub- 
lished many times in the open literature 1,3,6 as they 
have evolved. The flow is presumed to be governed by 
the unsteady Reynolds- aver aged Navier-Stokes equa- 
tions which are written in a body-fitted coordinate sys- 
tem. They are written in a usual conservation-law form 
as 

Q,t + {P — &v),t ? + {& ~ &v),{ = 0 

Here the subscripts with a comma denote partial differ- 
entiation, the subscript v identifies the viscous terms, 
and the superscript A indicates scaling with respect to 


the Jacobian of the coordinate transformation. Details 
for these terms are included in Appendix 1. With the 
ideal gas assumption, the pressure and total enthalpy 
can be expressed as 

p = (nr - 1- )(E-±pq 2 ) ; 

For the present study, the thin layer approximation to 
the governing equations is invoked 

& = = 0 

thus accounting for viscous flux terms only in the £ 
direction (normal to the body). 

Turbulence effects are accounted for through the 
notion of an eddy viscosity and eddy conductivity. 

li = fie + lit = Vi (1 + ihIih) 
k = kt + k t = kt (1 + k t /k t ) 

_ . IMPrA 

Pr e V ) 

The algebraic turbulence model developed by Baldwin 
and Lomax 12 is used to evaluate appropriate turbu- 
lence quantities. For separated flow regions, the no- 
tions of Degani and Schiff 13 are drawn upon to deter- 
mine proper turbulence length scales, A recent exten- 
sion to this fairly standard approach has been given 
by Hartwich 14 for massively separated and transitional 
flows. 

Computational Method 

Discretization of the governing equations results in 
a consistent approximation to the conservation laws in 
integral form 

where the time rate of change of the state vector Q 

within a cell is balanced by the net flux / across the cell 
surface. Flux quantities are represented using the up- 
wind biased, flux-difference-splitting approach of Roe 15 
with third-order accuracy. Salient aspects of this for- 
mulation are included in Appendix 2 and additional 
details have been given by Vatsa et al 6 . Solutions are 
advanced in time with a spatially-split, three-factor ap- 
proximate factorization method in diagonalized form. 

Interface quantities between longitudinal blocks 
are determined in a conservative second-order accurate 
manner. At the interface between the two blocks, the 
conserved variables are interpolated across the over- 
lapping mesh area using a technique introduced by 
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Ramshaw 16 . The conservation of mass, momentum, 
and energy are maintained across the interface bound- 
ary by redistributing the fluxes from one side onto the 
cell faces of the opposite side, A detailed discussion of 
the patching algorithm is given by Thomas, et al 17 . 

Surface Definition and Grid 

The surface definition for the F/A-18 forebody- 
LEX geometry was obtained from a detailed CAD de- 
scription of the 6% scale model tested by Erickson 10 , 
The configuration was modeled from the nose longitu- 
dinally back to $/ 5 = 2,337, a station slightly ahead 
of the LEX-wing juncture [x/d = 2,482) at which 
point the LEX leading-edge has already achieved 90° 
of sweep, Fig, 1. A downstream extension of the sur- 
face is then obtained by repeating the cross-sectional 
geometry from this station as described in the follow- 
ing paragraphs. 

The grid is divided into two longitudinal blocks 
which interface along a plane x = const* at the LEX 
apex station, Fig, 2, The forward block extends up- 
stream of the nose by approximately 1,5 5 and the aft 
block extends downstream of the LEX-wing station by 
approximately 1.5 5; both blocks extend radially from 
the model centerline by approximately 2,4 (For the 
downstream extension, the model cross-section is held 
constant.) The reference length er is comparable to the 
longitudinal extent of the forebody and 3/4 of the lon- 
gitudinal extent of the LEX-body juncture, 

The three-dimensional grid for each block is con- 
structed from two-dimensional O-type cross flow grids 
which are longitudinally stacked, constituting an H-O 
topology. Both blocks have 65 points in the radial di- 
rection; the upstream block has 31 points circumferen- 
tially whereas the downstream block has an additional 
34 points in the circumferential direction to represent 
the LEX, A more detailed view of the discrete sur- 
face representation is shown in Fig. 3(a). With the 
blocked approach, the LEX apex is represented as a 
singular point while the body representation is contin- 
uous through this station. Longitudinally, the grid is 
clustered near the nose and LEX-apex regions with 17 
stations on the forebody and 22 stations on the LEX- 
body. The grid is completed with 10 stations extending 
upstream from the nose and 9 stations extending down- 
stream from the station x/ff = 2.337; both extensions 
incorporate longitudinal stretching to provide resolu- 
tion near the configuration. (The grid upstream of the 
nose was treated as a separate block for turbulent flow 
computations). The F/A-18 forebody-LEX surface ge- 
ometry is represented with a total of 1,957 points, and 
the entire grid consists of approximately 185,000 points. 

The crossflow grids are generated using established 
transfinite interpolation techniques 18 * 19 with a method 


applicable to slender shapes 20 . A baseline grid (of ap- 
proximately 185,000 points) is first established with suf- 
ficient normal clustering near the surface to adequately 
resolve the laminar sublayer for the turbulent bound- 
ary layer flow at the subject wind-tunnel frees tream 
conditions (M^ = 0,6, J? e = 0,8 X 10 6 , and a — 20°), 
This grid produced an average normal cell she next to 
the wall of approximately 10 which corresponds to 
« 2 for the turbulent computations; a laminar sub- 
layer generally extends out to « 8,5. The same grid 
is used for the laminar computations as well. Fig, 3(b) 
illustrates the nearfield grid structure about a typical 
F/A-18 LEX-body cross section. This figure illustrates 
the challenge that is associated with generating a sin- 
gle O-type grid around such complex cross section with 
various break points in the surface geometry. Addi- 
tional details of the grid resolution in the LEX-body 
juncture as well as the LEX leading-edge regions are 
shown in Fig, 3(c) and Fig. 3(d), respectively. 

At the subject flight freestream conditions (Moo — 
0.34, R c — 13,5 X IQ 6 , and a — 19°) the baseline grid 
produced 8, Although an average t/ + slightly 

less than 10 has been reported 8 to provide adequate 
viscous resolution, a grid sensitivity investigation was 
conducted to provide solutions at flight conditions with 
comparable viscous resolution to those already achieved 
at wind-tunnel conditions. The baseline grid was mod- 
ified to have twice the number of points in the normal 
direction with increased clustering such that approxi- 
mately 4,5 fine grid cells were packed in the first cell 
of the baseline grid. The longitudinal and circumferen- 
tial grid distributions were not altered. The resultant 
fine grid consists of approximately 370,000 points and 
yields w 3 for the turbulent computations at flight 
conditions. 

Results and Discussion 

All computations were performed on the NAS Cray 
2, located at NASA-Ames. On this machine, the al- 
gorithm requires approximately 22 seconds per grid 
point per cycle. Converged results were obtained in 
nominally 2200 cycles requiring about 2.5 hours of com- 
puter time for the baseline grid. This number of cycles 
was sufficient to reduce the residuals by two to three or- 
ders of magnitude and to reduce oscillations in Cl and 
Cd to a negligible level. The computations were per- 
formed without the use of mesh sequencing or multigrid 
iteration. 

Results are presented for the subject wind-tunnel 
and flight test conditions. All comparisons of computed 
results (e.g, laminar vs. turbulent, baseline vs. refined 
grid, etc.) within a figure are presented from an iden- 
tical vantage point. The magnitudes associated with 
contour quantities are displayed with a color bar. 



Wind-Tunnel Computations 

The computed total pressure contours and surface 
streamline patterns are shown in Pig. 4 for both lami- 
nar as well as turbulent flow calculations. These solu- 
tions are obtained at M w — 0.6, = 0.8 X 10 6 , and 

a — 20°. The cross flow total pressure contours for the 
laminar computations (Pig, 4(a)) clearly show the LEX 
primary and the secondary vortex systems. In addition, 
the same crossflow contours show an appreciable body 
type separation under the LEX which will be discussed 
later in detail. Furthermore, there exists a clear pri- 
mary separation line (as indicated by the converging 
surface streamlines) on the forebody along with a sec- 
ondary separation line on the leeward side of the aft 
forebody. These are fairly flat bubble-type separations 
which, in general, occur on the order of a boundary 
layer thickness. The primary and the secondary sepa- 
ration lines merge and wrap around the LEX apex with 
a subsequent entrainment into the juncture between 
the LEX upper surface and the body. The turbulent 
flow computations demonstrate a drastically different 
flow structure as shown in Fig, 4(b), For example, it 
appears that the primary and the secondary forebody 
separations are completely eliminated in the turbulent 
case, Additional differences are discussed subsequently. 

A closeup view of the LEX-body flow structure is 
shown in Fig, 5 for both laminar and turbulent flow. 
These results clearly show the very complex flow that 
exists under the LEX, especially for the laminar flow so- 
lution (Fig, 5(a)), It appears that for this result, there 
is a bubble type separation under the LEX apex, sim- 
ilar to what occurred on the forebody. Further down- 
stream a larger vortex-like separation is evident along 
with secondary separation effects. Once again there 
are significant differences between laminar and turbu- 
lent computations as shown in Fig. 5(b). The body 
separation under the LEX is changed considerably by 
occurring closer to the LEX lower surface and practi- 
cally filling up the entire corner region. The separation 
pattern for the turbulent case is much .simpler than the 
corresponding pattern for the laminar case. 

Fig. 5 also provides a close view of the surface flow 
in the vicinity of the longitudinal grid-patching station 
(i.e., the LEX apex station). Note that the flow prop- 
erties are very smooth across this interfacing patch be- 
tween the two blocks of grid. 

The LEX upper surface flow pattern computed at 
the wind tunnel flow conditions for both laminar and 
turbulent flow are shown in Fig. 6. The laminar com- 
putation (Fig, 6(a)) clearly shows the secondary sepa- 
ration line with the subsequent reattachment line (in- 
dicated by the diverging surface streamlines) between 
the primary separation line and the leading edge. A 


tertiary separation is also indicated outboard of the 
secondary separation line. The turbulent computation 
(Fig. 6(b)) show the outboard movement of the sec- 
ondary separation line, as expected, and reduced evi- 
dence of the tertiary flow. 

The normalized density contours at various LEX 
crossflow planes are shown in Fig. 7 for laminar and 
turbulent flow. The color bar in the upper right cor- 
ner of the figure indicate the level of compressibility 
associated with the core of the LEX vortex system. In 
the vortex core region there is a 40% to 50% expan- 
sion compared to the freestream condition. Note that 
the laminar and turbulent computations produced very 
similar density contours. Additional analysis indicated 
that the majority of the core flow achieved a local Mach 
number of 0.9 or greater with a small supersonic zone 
occurring near the LEX apex. 

Contours of the static surface pressure coefficient 
are shown in Fig. 8 for both laminar and turbulent 
flow computations. The results show a compression re- 
gion around the apex of the LEX as well as the foot 
print of low pressure associated with the LEX primary 
and secondary vortices on the LEX surface. Although 
these surface contours of the laminar and turbulent so- 
lutions look very similar, a more detailed assessment 
indicates significant differences as will be discussed sub- 
sequently in conjunction with experimental wind tun- 
nel data. The experimental stations are highlighted in 
white and are located at xjz = 0.334, 0,587, 0.891, 
1,390, 1.701, and 2.143. For reference, the LEX apex is 
located at x/z — 1.006 and the LEX-wing juncture at 
x/c = 2.482. 

The comparison of the computed surface pressure 
coefficients for both laminar and turbulent flow are 
shown in Fig, 9 at the six different stations. The ex- 
perimental data are also shown for both the starboard 
and port sides of the model to assess flow symmetry. 
Figures 9{a)~(c) show the variation of the surface pres- 
sure coefficients as a function of azimuthal angle 6 on 
the forebody. The windward and leeward sides of the 
forebody correspond to 9 = 0° and 180°, respectively. 
Both laminar and turbulent solutions slightly under- 
predict the pressure peak at station 1 as shown in 
Fig. 9(a). This underprediction of the pressure peak 
could well be attributed to the inability of the H-0 grid 
topology to resolve the flow near forebody nose. This 
figure also indicates the differences between the lami- 
nar and turbulent solutions which occur in the range 
130° < 8 < 165°. This difference is due to the sep- 
arated zones on the forebody which occur for laminar 
flow. Figures 9(b)-(c) show a very good correlation 
between the computational results and the experimen- 
tal data* particularly the turbulent flow solution. The 
differences between theory and experiment shown in 



Fig. 9(c) near 6 = 0° and & = 180° are believed to be 
associated with circumferential grid resolution. 

The LEX computed surface pressure coefficients 
are shown in Fig. 9(d)- (f) along with the experimen- 
tal data as a function of LEX exposed semispan, yfs . 
In this way, yja — 0 corresponds to the LEX-body 
juncture and yjs = 1 corresponds to the LEX leading 
edge. Figure 9(d) reveals a good agreement between 
the experimental data and the turbulent flow computa- 
tion. This figure also shows that the laminar flow under 
predicts the low pressure peak that is associated with 
the primary vortex system. This difference between 
the laminar and turbulent solutions was not as evident 
on the color surface pressure contours shown earlier in 
Fig. 8. The correlation between theory and experiment 
begins to degrade at station 5, Fig. 9(e), and continues 
to do so at station 6, Fig, 9(f), It is believed that this 
difference is primarily associated with wing-induced up- 
wash effects present in the experimental data. The 
test was conducted with a complete F/A-18 configu- 
ration, including flow-through inlets, whereas the com- 
putations only modeled the forebody-LEX portion of 
the configuration. The wing- induced effects would be 
strongest near the aft portion of the LEX and, in gen- 
eral, would cause more negative suction pressures by 
strengthening the LEX vortex. Inlet-induced effects are 
a secondary source for this discrepancy. None the less, 
it is noteworthy that the forebody pressures as well as 
the pressures on the forward portion of the LEX can be 
accurately predicted with the turbulent computations 
on the isolated forebody-LEX geometry. 

The correlation between theory and experiment is 
further assessed by examining the experimental sensi- 
tivity of the surface pressures to changes in angle of 
attack as well as Mach number. Fig, 10. In this fig- 
ure results are presented for stations 3 and 4 which 
represent experimental perturbations about the subject 
wind-tunnel conditions; the turbulent theory from the 
prior figure is also included. The data are shown to be 
quite sensitive to changes in angle of attack of ±2.5° as 
well as to changes in freestream Mach number from 0.4 
to 0.7. The differences between theory and experiment 
are, in general, significantly less than the the incremen- 
tal effects of a and jWoo shown in the experiment. 

Flight- Test Computations 

The flight-test solutions are obtained at = 

0.34, = 13.5 X 10 6 , and a = 19°, Due to the high 

Reynolds number associated with flight, grid sensitiv- 
ity effects were assessed. Various flow parameters (e.g-, 
Mach contours, density contours, surface flow patterns, 
surface pressure contours, etc.) were contrasted be- 
tween the solutions obtained with either the baseline 


or the refined grid. No major differences between these 
solutions were found. Typical results which contrast 
the baseline and refined grid solutions are presented 
in Fig. 11 for static surface pressure contours and in 
Fig. 12 for the LEX upper-surface flow pattern. It was 
concluded from this study that the baseline grid pro- 
vided credible solutions which could be used for the 
subsequent correlations with the flight test results. 

The total pressure contours at various crossflow 
stations as well as the surface flow pattern computed 
at flight condition are shown in Fig. 13. Fig. 13(a) 
shows an overall view of the flow, whereas Fig. 13(b) 
illustrates the LEX-body flow details. The fore body 
surface flow pattern qualitatively resembles the turbu- 
lent solution obtained at wind tunnel flow conditions 
(c.f., Fig. 4(b)) in that the flow appears to remain at- 
tached. The LEX-body flow exhibits greater differences 
between flight (Fig. 13(b)) and wind-tunnel (Fig. 5(b)) 
conditions principally with respect to the body sepa- 
ration under the LEX. This separated zone is smaller 
at the flight conditions than was observed in the wind- 
tunnel computations. 

A side view of the computed surface flow stream- 
lines on the forebody at flight condition is shown in 
Fig. 14(a). This figure also shows part of the LEX ge- 
ometry. The surface flow pattern obtained in flight is 
shown in Fig, 14(b) from a similar vantage point. The 
body separation line under the LEX is well predicted 
by the theory. Furthermore, the overall surface flow 
pattern around the forebody agree well, qualitatively, 
with the flight test results. Similarly, Fig, 15 show a 
correlation between the computed surface streamlines 
and the flight test results on the windward side (of the 
front part) of the forebody. The correlation between 
these results is also quite good. 


Concluding Remarks 

Three dimensional viscous flow computations for 
the F/A-18 forebody-LEX are. presented for both wind 
tunnel as well as flight conditions. At wind tunnel con- 
dition, significant difference between laminar and tur- 
bulent solutions are revealed, particularly on the sur- 
face flow pattern. Good agreement between the com- 
puted surface pressure coefficient and the experimental 
data are also disclosed. The turbulent computations 
provided a better correlation with the data. At flight 
condition, the grid sensitivity study revealed that the 
baseline grid clustering with only 185,000 points was 
sufficient to adequately resolve the viscous laminar- 
sublayer for the turbulent computations. Furthermore, 
the computed surface flow pattern correlated well, qual- 
itatively, with the flight test results. 
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Appendix I - Governing Equations 
The in viscid flux terms from the governing equa- 
tion are defined as follows; 
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The contravariant velocity components in these 
equations are given by: 

+ (, V V + 

The viscous flux terms from the governing equation 
are defined as follows: 
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where 



u. ^Ci + Z,x<f>t 2 

v, t4>t x + £,v<t>t 2 

l ((%•).« + Pr^-l))^ei + U ^t2 J 

0 

v ,n$n l + 

i ((V).n + J 
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u ,f^fl + ?,V^{2 

l ((%•).? + +^^f2 J 


(*•)* + (O) 2 + (U* 
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The transformation metrics are: 


f , x V,z t y x 

£.v WtV ( t y 

*1,2 5,2 J 


- *,nV* *AV* - V,t*« &<*.* - *,«&, 

*,**,* - “ *,«** “ *,e*,n 
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and the Jacobian of the transformation is given by: 

J - 1 « *,<(&„*,* - z,„y,f) - y,<(*,t,*,f - * t ** tf ) 

+ *,«(*■*».* 

Appendix 2 - Upwind Algorithm 

Spatial derivatives are expressed as a conservative 
flux balance across a cell such as 

(Sf/S()i = (V* - 4- ./,)/(%+ 1* - fi-1*) 

and the interface flux is expressed as 

= V2[^(Ql) + /(9ji)- I i I (Q« “ «L)li+i A 

Here Q/, and <3# are the state variables in the left and 
right cells and 

A s ar/SQ = T(A + + A -Jr -1 

T and T -1 are diagonalizing matrices and A is the di- 
agonal matrix of eigenvalues with 

A* = y 2 (A± I A I) ; I A |= T | A 1 3T 1 

The superscript denotes evaluation of the matrices 
with Roe- averaged quantities such as 


so that 


P = \[PLPR 

U = [ui + U R \/p R / pl )/{1 + « 

5 = + «nVpn/pL)/( 1 + Vpr/pl) 

t® — (u>z, + wjjVPb7ac)/(1 + Vpr/Pl) 
6 = (/fi, + i¥ rs/pr[pl ), / (1 + V Pr/pl) 
a 2 b ( 7 - l)(# - (ft 2 + 0 3 + ui 2 )/2) 


is satisfied exactly. 

State variables at the cell interfaces are formed 
from interpolation of the primitive variables g = 
[p, u, v, tu, p] T by the one-parameter family 

(®0 <+ i£ = « + V«[(i - K )v?» + U + 

(««),-+ i/j = ?i+i - V4(l + «)V9< + i + (1 - «)Ag<+i] 
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where the backward and forward differences are ex- 
pressed respectively with a typical minmod limiter as 

y g — mmrnod[y g, BAg] 

A q — minmod[Aq , B y g] 

with 

minmod{x, y) ^ sign{x)x 

max [o,min(a; stgnfy), y stgn(a:))] 

Here B is a compression parameter 

B = (3 — «)/(l — «) 

and k — 2 /3 for third-order accuracy. 
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(bj Surface representation. 


Figure L- F/A-18 HARV vehicle. 
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(b) Ns a? fie id* 

Figure 2,- Blocking strategy for F/A-1S forebody-LSX grid, 
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(a) Laminar flow. 



Figure 4,- Total pressure contours with surface Sow pattern - overall view. =0.6, K s = G.S X lG a ,a = 20° 


(a) Laminar flow. 




(bj Turbulent Sow. 

Fi^uro sj ■■ Total pressor s? contours with surface Sow pattern - Lflf.X-body detail. ■&£■>:; — 0,Xi <A * 0>B v* ItJ^os — - *p0 > 
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(a) Laminar flow. 



Figure 6,- LEX upper-surface flow pattern. Moo — 0 * 6 , J? c — 0.8 x IQ 6 , a = 20 °, 
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(a) Laminar flow. 




(b) Turbulent flow. 

Figure 8.- Static surface pressure coefficient contours. M = 0.6, — 0.8 X 10 6 , a — 20 Q . 
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Figure 9.- Correlation between theory and experiment - surface C p . ™ 0.-6, R c = 0.8 x 10% a =■ 20°.. 
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(a) Station 3, Moo = 0*6. 




(c) Station 4, M OT = 0,6. 



(d) Station 4, a ™ 20°, 


Figure 10.- Experimented sensitivity of surface Cp. J?* = 0.8 X 10 6 . 
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(a) Baseline grid* 



(b) Refined grid. 

Figure 12.- LEX upper-surface Sow pattern. — 0.34, J? e — 13.5 x XO r> ? a ™ 19^. 
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[a) Overall view. 



|b) U£X- foody detail, 

igure IS.- Total pressure contours with surface Bow pattern. ----- G..M, E c - IS, 5 x a .19* . 







(a) Computational result. 



(b) Flight test, 

Figure 15.=- Comparison of computed turbulent surface streamlines with flight test - bottom vie tv. 

Afoo - 0,34, = 13.5 X 10®, a «* 10°. 



